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The two-dimensional generalization of the Polyelectolite model of DNA is proposed. 
■ It is reduced to the boundary problem for nonlinear quite integrable equation sinh- 

Gordon. In the linearisible version the exact solution is constructed and its asymp- 
totic is found. The soliton solutions of nonlinear equation are calculated that allowed 
tells about the possibility of the structural phase transition in the considered system 
(DNA+polyelectrolite) on the temperature. 
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At the present time constructing and investigating mathematical models of DNA is 
an important problem of both theoretical biophysics and mathematical physics. Modern 
achievements in the field of nonlinear equations theory that occur, as a rule, in these 
models allow not only to set initial and boundary problems for such equations, but also 
to solve them in a number of cases. 

The given work is devoted to construction of exact solutions for a two-dimensional 
generalization of so-called polyelectrolyte model of DNA (PM DNA) [1]. In essence here 
C\| ' we deal with solution of an inverse problem in the broad sense, i.e. it is to judge structure 

and properties of a DNA molecule from a field restoring within the framework of this 
l/-) , model and also values of some functionals of it. 

It is necessary to note that PM represents the most concentrated description of role 
of Coulomb interactions which are of a great importance in formation of structure of a 
DNA molecule and in its functioning as well as in intrinsical structural transitions. One 
more important advantage of the model is possibility to operate with terms of statistical 
physics, which allows to implement a clear mathematical problem statement. 

The formulation of model and statement of the problem. In reality a DNA 
molecule is a strongly charged poly-ion, with two charges of electrons falling at every pair 
of nucleotides. The idea of PM consists in replacement of such system by an infinite and 
regular negatively charged cylinder of radius vq with a given surface charge density placed 
into solution of polyelectrolite. It is supposed that a processes of relaxation proceeds 
quickly enough, so only an equilibrium situation is considered. Besides here we are re- 
stricted ourselves by consideration of the case when a solution consists of unitary charged 
ions and electrons. 

Standard physical reasons result to a system of Poisson-Boltsmann equations written 
down at any point of volume outside the cylinder: 

div E = 4np, E = —grad U, p = e(rij — n e ), (1) 

in this case n ifi being under condition of thermodynamic balance are defined by Bolts- 
mann's distributions: 
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1 The work is executed at financial support of the Russian Found of Basic Researches (project No 
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ni(U)=n 0i e r ? n e (U) = n 0e e t . (2) 

In formulas (1), (2) the following notations are used: E = E(x, y,z) is a vector of elec- 
trical field intensity, U = U(x, y, z) is a dimensional scalar potential, e is a charge of 
electron, p is a charge density, rii, n e are concentrations of ions and electrons at the 
point r = {x,y, z} correspondingly, n 0i , n 0e are equilibrium concentrations, T is absolute 
temperature, Boltsmann's constant being assumed to be equal to unit. 

Supposing that n 0i = n 0e = n (the condition of total system neutrality) as well as 
that there is a symmetry in shifting along an axis z being a symmetry axis of the cylinder, 
we will pass the plane orthogonal to this axis. Then introducing dimensionless variables 
r' = y/(2tb)i u = e U/T, where is Debye radius, r 2 D = T/(8iren ), from (1), (2) we will 
obtain a basic equation of the model (strokes at independent variables are omitted) 

Au = 4 sinhw, (3) 

were A is two-dimensional Laplace's operator. 

It should be noticed, that the equation (3) also can be obtained from a strict and 
valid procedure of breakage and disengagement of BBGKI chain in neglecting fluctuations 
[2,3], with the function u being an average potential of the system which is a functional of 
Coulomb's two-partial potential. The equation (3), as applied to PM DNA, corresponds 
to the case of a salt excess in a polyelectrolyte solution []. 

sinh -Gordon equation (3) belongs to the class of two-dimensional completely inte- 
grable equations of the elliptical type, for which representation of zero curvature is known 
[5, 6] and, hence, the Inverse Scattering Transform Method (1ST) can be applied. In the 
papers [7,8] this equation was solved for a half-plane y > 0, in [9] an 1ST procedure was 
employed in solving sinh-Gordon-equation fitting Coulomb's system on a plane at nega- 
tive temperatures, and in a series of papers [10-12] an equation of the same type (with a 
function sin ii) was solved on a whole plane with given conditions on chosen beams (for 
the purpose of obtaining unambiguous solutions). 

Further it will be convenient to work with the equation (3) rewritten in polar coordina- 
tes: 

1 1 

u rr + -u r + — u w = 4 sinh u, (4) 

where x = r cos (p, y = r sin tp, r > tq (here tq is a dimensionless radius of the cylinder) 
0. Let us assume that u — u(r, <p>) is a real- valued function, u —>■ with r — > oo quickly 
enough, with u(r, <p) = u(r, <p + 2n). 

For (4) are natural two statements of boundary problems: 

a) Dirichlet's problem for an exterior of circle: 

u\ r=ro = f x (ip), (5) 

2 If the right part of (3) to substitute for an exponent (having a negative sign both before its argument 
and before the exponent) will have Liouvillc's equation having an exact solution [4] that corresponds to 
a salt lack). 

3 In the one-dimensional case recently this equation has been numerically solved in [13]. 
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b) Neumann's problem for an exterior of circle: 

U r \ r =r = h{^>), (6) 

in which fi((p), hiv) are given periodic functions (linear densities of charges of the cylinder 
revealing a specific character of DNA molecule). 

Linearized model (a two-dimensional generalization of Debay-Chukkel's 
model). At r >> r the equation (4) can be linearized and then can be written as 



i 1 

7J" u>(f(p 



r r 

Taking u = u(r, ip) = R(r)&((p), we shall have 

r 2 R rr + rR r - (4r 2 + X)R = 0, $ w + A$ = 0, (7) 
where a parameter A G C. Considering a finiteness of solution from (7) we find: 

u(r,<p) = J2b m K m (2ry m ^ A = m 2 . (8) 

m 

where b m are constants, K m (.) are modified Bessel functions (McDonald's functions), and 
the summation is over m within limits from —00 to +00. In the case of the problem (5) 
coefficients b m are defined from the condition 



Thus, we obtain 

Similarly in the case of a problem (6) we find 

u(r, ¥) = J2 c m K m {2r)e m \ A = m 2 , 

m 

and the coefficients c m are defined by relations 

1 f o 27r dkpMtfe-*"* 
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27r^ m _ 1 (2r ) + ^ m+1 (2r )' 
Hence it follows that 



It is easy to show that both (9) and (10) meets the requirement of reality Representa- 
tions (9) and (10) are solutions of Debay-Chukkel's two-dimensional theory corresponding 
to anisotropic medium (polyelectrolyte) , which can arise owing to, for example, possible 
(probable) fluctuations in the system. 
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Fourier- factors fx m , / 2m of expansion of functions f\{f), f2{ ( fi), are obviously connected 
by the relations: fi m /K m (2r ) = f2 m /{.'2K' m (2r )), which allows to establish a nonlocal 
relation Q): 



u r (r , if) 



2~ 



df'H(r , ip - f')u(r , f) 



'ID 



in which the function H(r , if) is determined by the equality: 



7T ' " 



K' m {2r) 



71 m K ™( 2r o) 

Using (9) and the Fourier-factors connection (11) established above we obtain 



u(r,<f) 



An 



2ir 

dip'J2 eim(iP ~ V ' )K m(2r) 



h{iff) + h{f') 



K m (2r ) 2K' m (2r ) 



Returning to (8), where 



b m = — 
2vr 



2- 



dipe~ imv 

and considering that 6„ m = b m , we have 



u(r ,if) u r (r ,if) 



K m (2r ) 2K' m (2r ) 



u(r, if) = u {r) + 4 ^(Re b m )K m (2r) cos 



mip, 



;i2) 



m=l 



with «o(r) = 6o-^o(2r), bo = b . The formula (12) is convenient for obtaining asymptotic 
representations. Then 



u Q {r) 



'0 -2r 

e , r — > oo. 



Thus, we can deduce a relation describing an effect of Debye screening. Putting 

u + {r,f) = Y,b m K m {2r)e m ^ ~b m = 4Re b m , 



m=l 



and assuming the possibility of analytic continuation of b m on an index m, we convert a 
last sum in integral in the ordinary way: 



u + (r,f) 



ds ^- K s (2r), 



(13) 



l c 2isin7rs 

where a contour of integration starts from an infinite (along Re s) point located in the first 
quarter of the complex variable s plane, then goes parallel to the real axis crossing it at 
the point < Re s < 1, and then ends at infinity parallel to the real axis and remaining 
in the fourth quarter. 



The idea of the approach belongs to V.D.Lipovsky. 
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To calculation an asymptotic form of the expression (13) it should be applied the 
saddle-point method. For that at first we should notice that the function K s satisfies the 
equation 

K s , rr + -K a , r - (4 + — )K S = 0. (14) 
Taking £ = 2r, K s = f s /V£, we reduce (14) to the form 

2 1 

m — — 

f'J - (—^ + 1)/- = o. 

From here it follows that 



3~T 

' ~4 



or, after simple calculations (assuming s/£ = o(l), £ — > oo), 

1 s 2 

K s (2r) ~ e~ 2r e~^, r — > oo. 



Then it is possible to rewrite (13) as: 

M+(r )V 9) ~ ^e- 2r / ds — - e "( w +v)-^, r^oo. (15) 

v 7 ^ 7 7c 22sinvrs ' v ; 

Let us assume that s = so£, so = 2i(7r + </?)r. Then a whole picture will be turned 
through a corner - ir/2: the poles will turned out in the negative part of the imaginary axis, 
and the integration contour C will turn into a contour C beginning at the point a — too, 
going parallel to the imaginary axis and then crossing it and finally ending at the point 
—a — ioo remaining parallel of the axis, where a is a positive constant. Assuming that 
"factors" b s permit an analytic continuation to required areas of the lower semiplane, 
instead of (15) we obtain 

"2r e „ f ^ 6 ( s oQ r -r(v+ V ?(2j-e) 



On introducing a function = £ 2 — 2£ we find the saddle point: £ = 1; then in its 
vicinity = — 1 + (£ — £ ) 2 + — £) 3 ). Deforming the contour in such a way that 
it goes through the saddle point in the direction of the steepest descent, we have 

u+M = fJWLpW (i + o(I)) j e -~ir, r _ oo. 

I smh[2r7r(7T + </7jJ r J 

This expression, which is the leading member of the asymptotic form, gives an amend- 
ment to the effect of Debye screening on plane. 

Exact solutions and "soliton" configurations of the equation (4). Let us 

proceed to constructing exact solutions of the equation (4). It should be notice that the 
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solution of the boundary problems (4), (5) or (4), (6) encounters serious mathematical 
difficulties involving application of 1ST (a scattering problem for the operator of an as- 
sociated linear problem on half-axis). Therefore here we should be restricted by a more 
simple problem: in neglecting the radius tq we will construct exact solutions of (4) on a 
whole plane by means of Darboux Transformation (DT) method. 

Employing a direct calculation it is possible to check up that the equation (4) is a 
condition of compatibility of the overdetermined linear matrix system 

^ r = U^, V v = rW, (16) 

where * = *(r, (p, A), U = U(r, ip, A), V = V(r, ip, A) G Mat(2, C), A G C is a parameter, 
and the matrixes U, V are given as 0: 

U(r, „, A) = -e-<x 3 + -^e ^—a 2 - — e ^, 

y(r, *>, A) = - ^e"'V 3 - - 

Here <7j, i — 1,2, 3, are Pauli's standard matrixes. 

Further instead of variable r, 9? it is convenient to introduce the "conic" variables ( = 
(ip + iha.r)/2, ( = {ip — i\nr)/2, and also using invariance of a condition of compatibility of 
the system (16) relative to a cyclic rearrangement of Pauli's matrixes, it is worth changing 
to the another its gauge. Then instead of (16) we have the following 2x2 matric system 

^ c = AV, * f = BV, (17) 

where 



The compatibility condition (17) becomes Aq — + [A, B] = 0, and an appropriate 
nonlinear equation, takes the form 

= 4 sinh u e ~ 2l(c ~^ , 

which, as it is easy to show, is equivalent to (4). 

Let \& = (vl/t 1 ), \I/( 2 )), vI/W = (6,x) T - Then the equations system (17) can be written 
as a system of four equations 



Let 9%, Xi be a fixed solution of (18) corresponding to the choice A = Ai. Let us assume 
5 For the first time these expressions were obtained by S.S. Nikulichev 
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t>i Xi 

and check up covariance of the system (18) relative to DT of a such form. After simple 
calculations it can be obtained 

M = M + 21n^i. (19) 

In the theory of DT the relation is called a "dressing" one. 

To construction an explicit solution it is necessary to assign some initial one; here as 
the solution we choose u — 0. Turning back to polar coordinates, then from system (18) 
we obtain 

6i r = -^t[cos tp(^- - Ai) - % sin + Ai)]xi, 
li Ai Ai 



(20) 



The given system has an integral: Of — Xi = Af, where Ai is an arbitrary, generally 
speaking, complex constant, which allows to integrate (20). We have 

Xl = l Al {e T ^ hlM - e- ri+lnAl ), 6 1 = ^i(e ri - lnAl + e" ri+lnAl ). 

Here r x = —1/(2?) [cos </?(l/Ai — Ai) — isin<£>(l/Ai + Ai)]r + ln_B 1? Bi is a real constant ( 
for simplicity A 1 is considered to be real as well). 

Considering (19), it is can be seen that the potential becomes real at Ai = e iai , cti £ 
R, cti G [0, 2ir]. Then the simplest solution of the equations (4) takes the form (u[l] = u) 

u[l\ = 21n{tanh[rsin(v? + ai) + (21) 

where 5\ = In The solution (21) is an analogue of the solution type of 1-kink for 

nonlinear equations of the hyperbolic type [5] and describes a distribution of the potential 
on plane. Also it is of interest to note that in solving the equation (3) on half-plane by 
means of 1ST, 1-kink solutions turn out to be forbidden because of the potential must be 
real [8]. 

Let us consider an double dressing procedure for the initial solution. Repeating pro- 
cedure of obtainment (19) we have 

" |2| = * + 21n W (22) 



Here 



9 2 [1] = e - 2 <"(0 2C - - -£o 2 ) = X2X2 - Ax^X 

X2 [l] = e- 2 *(x2< - ^2) = \2h ~ \1-X2, 
Xi Xi 
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and X2 = X2( r , #2 = d 2 (r,<f) are solutions of the system (20) for A = A 2 . Then the 
relation (22) can be rewritten in the form 

roi 1 01 ^2X1^2 — Ai#ix 2 / oq n 

* l= " + *w^S' (23) 

On examining realty of the expression (23) it can be resulted in two possible cases. 

1. Let Afc = e lc * k , a k G [0, 27r], k = 1,2, i.e. both complex parameters must belong 
to an unit circle 0. Let us assume Tk = r sin(<£> + a k ) + ln_Bfc, B k G R, k — 1, 2, and as 
solutions of the system (34) take 

Xk = C k e™x + iD k e- Vk - l * & \ 9 k = C k e™ k - iD k e- T ^ 5 K 

In this case the realty is guaranteed provided Xk = 9 k , C k , D k E R, —4iC k D k = A\, 
and then we obtain: 

m o^ l + hks h h( \ sinh(r 2 - Tx + z/x) 

u ks [2 = 21n i r— > h k3 = h(r,<p) = -cot Q — -. (24) 

1 - /i fcs cosh(ri + T 2 + v 2 ) 

Here ^ = ln|5 2 /5i| - (1/2) In l^i/i^l, z/ 2 = In |<5i5 2 | + (1/2) In l^^l, 5 1; <5 2 G R, i4 = 
A/Cfc, q=( ai -a 2 )/2. 

The solution (24) is an analogue of 2-kinks solution for equations of the hyperbolic 
type. It differs a little from a solution of the equation (3) obtained in [8] for half-plane 
by means of 1ST (ibidem the appropriate asymptotic forms at r — > 00 have been found). 

2. Let Ai = 7ie iai , A 2 = 1/A X , <*i, 71 G R, <*i G [0,2vr], 7l > 0. Then r(Ai) = 
f (A 2 ), r(Ai) =Vx = -l/(2i)[cos^(l/Ai - Ax) - zsin^l/A! + Ax)]r = T 1R + iT n , T 1R = 
(l/2)siii(^ + ai)(l/7i + 7 1 )r, r ir = (1/2) cos(</? + «i)(7i - l/li)r. Setting Q x = x 2 , 9 2 = 
Xi, k = 1, 2, where Xfc — F k e k — G k e~ Tk , 9 k = F k e Tk + G k E~ Tk , we have the other real 
solution: 



roi i\ l ~ hb u ut \ + n \ sin[2cos((^ + ai)(sinh7 1 )r] 

u b [2 = 21n — — , h b = h{r,ip) = cot (In 71) — , — ■ - — — -. 25 

1 + cosh[2 sm{if + «i)(cosh 7l )r + n x \ 

Here G fe G R, 4F k G k = A 2 k , A\ G R, m = (1/2) ln[(Ff + + 

The formula (25) is an analogue of the solution type of breather (a double soliton) 
for a hyperbolic equation. In [8] the solution has been obtained by help of 1ST and its 
asymptotic forms have been calculated. 

Procedure of obtaining subsequent "dressing" relations is same as in deducing (23) 
(see, for example, [14]), and also results in N - "soliton" solutions. 

On Phase Transition in PM DNA. In the context of the problem under examina- 
tion a rather wide solutions spectrum calculated above can be interpreted as one answering 
to various parameters or, perhaps, to various states of DNA. It allows, in particular, to 
put forward a hypothesis about a probable (within the framework of a used model) struc- 
tural phase of transition of the 2-nd sort on temperature (for example, type of spiral — > 
ball or B-form — > Z-form, where B and Z are right and left spiral DNA accordingly). Such 
a transition is possible due to our system is two-dimensional one (it should be reminded 

6 Employing 1ST for half-plane it has been shown that a parameter A has sense of a spectral parameter 
for an associated linear problem (16) [7,8] whereas values X k form its discrete spectrum. 
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that, according to Landau [15], such changes are forbidden in one-dimensional systems). 
As its mechanism the following one can be considered. 

Let two complex parameters (an eigen values) forming an inversion relative to an 
unit circle in the plane of a complex parameter A (state of the "breather" type) are 
available. It is assumed that in the system some change of temperature occurs, so that 
given parameters become dependent on it. It means that they start moving along the 
beam drawing between them (on the complex plane). In bringing together in the vicinity 
of the unit circle "a scattering" of the parameters occurs one on the other, and as a 
result they either with some probability p = p(a,T) drift apart along the circle, or with 
the probability 1— p "pass" through each other in such a way as to keep forming inversion 
relative to the circle obeying requirement of realty of the potential. The temperature, at 
which the phenomenon take places, is critical for the process; and the process corresponds 
to a structural phase of transition of the 2-nd sort Q In order to prove this, following the 
idea of Landau, we are to calculate free energies of 2-kinks and breather configurations 
in the vicinity of the critical temperature and obtain a (nonzero) jump of the second 
derivative with respect to a free energy on temperature. 

However, there is one essential difficulty The point is that according to so-called 
Derrick's theorem, in two-dimension systems a functional of energy turns out to be infinite 
(see [16], where questions on stability of integrable equations solutions have been detailed 
considered as well). Therefore to estimate the free energy functional, strictly speaking, 
it is to employ some mathematically correct procedure. Nevertheless, we suggest quite a 
simple (and correct, in our opinion) technique allowing to avoid this difficulty. 

Let us identify a free energy with a generating functional of the model (4): 

F(u) = [ [ ( 1 dxdyU(coshu-l)) + ^(ul + u 2 y ] 




= J j rfr^r{4(cosh M -l) + i( Mr 2 + !<)}, 

where a symbol jlR 2 \ C r=ro | means that one integrates over the space M. 2 with a remote 

circle of radius r$. Using the representations (24) or (25) for a function u = u(r,(f), we 
reduce J-'iu) to the form convenient for subsequent estimations: 



x Ah 2 + h 2 + ^ 



Here an auxiliary function h = h(r, <p) satisfies a nonlinear (and also, as well as (4)) a 
completely integrable equation 

h ■ 2 ^r + % , K 4fe(l + h 2 ) 

rr+ r 2 w + l-h? + r ~ l-h? ' 
whose representation of a zero curvature can be obtained from the conditions of system 
compatibility (16). 



7 In other words, here the phase of transition is connected to an asymptotical degeneration of eigen 
values, and the potential plays a part of the ordering parameter. 
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According to the critical phenomena theory [17,18], close to the critical point it is to 
assume 71 = 71 (t) = e Vl , where V\ = Vi(t) ~ t dl , t = (T* — T)/T*, di is the critical factor, 
di > 0, T* is the critical temperature. Then a simple estimation for h = h b near T — T* 
gives 

1 sin[2 cos(^ + ai)(7ir] _ 1 

~ ; = —F b (r, (p), v\ — > 0. 

vi cosh [2 cos(y? + a x )C 2 r + /ii] f 1 

Here Ci, C 2 is a real constant. 

Operating analogously for a condition characterizing by two "kinks", on setting v 2 = 
(«i — a 2 )/2 ~ t d2 , where d 2 a critical factor also, d 2 > 0, we find ( $ = y? + («i + «2)/2) 

1 sinh Ui 1 

/ife S ~ 7 : — -r-775 — ■ ■ j r = —F ks {r, <p), v 2 ^ 0. 

f 2 [cosh z/ 2 + smh(2r sin smn z/ 2 J f 2 

Thus, from (26) one can obtain 

Fb(t, r o) ~ vlQ b (r ), T ks {t, r ) ~ v%Q ks (r ), 

where Qb(r ), Qks( r o) are expressions that can be resulted in by what ever way of regu- 
larization of corresponding integrals. From here it follows that 

{.Fb)TT _ Qb di 2d x - 1 2 (<fa-<il) 



\J~ks)TT Qks d 2 2d 2 — 1 

Assuming d 2 — d\ = 0(t) at T /* and eliminating an uncertainty available here, we 



have 



{<Fb)TT\T=T* 7^ [3~ks)TT\T=T* ' 

and, thus, we show that at the critical point the existence of a change of phase of the 
2-nd sort is possible (this qualitative conclusion, certainly, does not depend on a way of 
regularization). It should be also noticed that exact values of the critical factors can be 
determined on obtaining identities of traces arising in analyzing an appropriate spectral 
problem (in [8] the values for half-plane have turned out to be equal to 1/2.). 

Conclusion. In the present paper an exact solution for Debay-Chukkel's two-dimensio- 
nal model has been constructed and its asymptotic form has been found as well. In the 
nonlinear case by means of Darboux Transformation method exact solutions have been 
calculated and a hypothesis about an possibility of a change of phase has been put forward 
and checked. Note that we did not solve a boundary problem for the equations (4), as 
such, because 1ST, generally speaking, is not adapted for this procedure Q. Therefore a 
conclusion about existence of a change of phase has, as a whole, a preliminary character. 
Besides a number of other important questions arising in the framework of PM such as 
question about thermodynamics of investigated system, about changes of phase corre- 
sponding to more complex soliton solutions, about methods of functionals regularization 
and others have not been considered either but will be discussed in future. 

8 Recently some progress in this direction has been achieved in [19] for a hyperbolic version of sin 
-Gordon equation. 
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